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ABSTRACT 

The  authors  present  the  extension  of  a previously 
developed  boundary  element  method  to  predict  the 
flow  inside  waterjet  pumps,  including  the  effects  of 
sheet  cavitation  on  the  blades.  The  circumferentially 
averaged  interaction  between  the  rotor  and  the  stator  is 
accounted  for  in  an  iterative  manner.  The  method  is 
applied  in  the  case  of  an  actual  waterjet  pump  and 
comparisons  of  the  predicted  and  the  measured  rotor 
torque  are  presented. 

INTRODUCTION 

Due  to  the  demand  for  high  speed  vessels,  the 
application  of  waterjet  propulsors  on  commercial  and 
navy  vessels  has  increased  in  recent  years.  Waterjets 
are  the  propulsion  of  choice  for  high-speed  naval  ships 
and  fast  ferries.  Compared  to  conventional  propellers, 
waterjet  propulsors  provide  several  advantages. 
Waterjet  propulsors  improve  maneuverability,  reduce 
the  possibility  of  cavity  occurrence  by  controlling  the 
flow  inside  the  casing,  and  reduce  the  likelihood  of 
damage  to  the  blades  by  protecting  them  inside  the  hull. 
Nevertheless,  like  in  many  other  fluid  machines,  such 
as  water  turbines,  pumps,  and  marine  propellers,  the 
performance  of  waterjet  propulsion  systems  is  affected 
by  cavitation  in  many  significant  ways.  The  foremost 
hydrodynamic  issue  is  the  thrust-breakdown  due  to 
cavitation.  This  inability  to  increase  thrust  is 
accompanied  by  noise,  vibration,  and  erosion.  The 
issue  is  more  serious  and  harder  to  address  in  the  case 
of  waterjet  propulsion  systems  than  in  conventional 
propellers,  inasmuch  as  waterjet  propulsion  systems 
have  several  components  such  as  intake  duct,  rotor, 
stator,  shaft/hub,  and  casing.  Due  to  the  interactions 
among  the  flows  around  the  various  components,  the 
flow  is  more  complex  and  multi-featured.  Presently, 
our  capability  to  predict  cavitation  and  its 


consequences  in  waterjet  propulsors  is  seriously 
limited.  Model  and  occasional  full-scale  measurements 
are  the  means  that  are  most  often  resorted  to. 

Numerical  methods  for  the  prediction  of  performance 
and  design  of  the  waterjet  rotor  and  stator  components 
were  presented  in  Taylor  et  al  (1998)  and  Kerwin  et  al 
(2006).  These  methods  were  based  on  inviscid  flow 
methods  (vortex-lattice  methods)  applied  on  the  blades 
of  the  rotor  or  stator,  coupled  with  either  Reynolds- 
Averaged  Navier-Stokes  (RANS)  or  Euler  equations 
solvers  for  the  solution  of  the  global  flow  through  the 
pump. 

Chun  et  al  (2002)  and  Brewton  et  al  (2006)  applied 
RANS  methods  on  the  rotor  and  stator  blades,  where 
the  interaction  between  the  rotor  and  the  stator  was 
considered  in  the  imsteady  sense  by  the  former,  and  in 
the  circumferentially  averaged  sense  by  the  latter. 

A comprehensive  review  of  issues  concerning  the 
prediction  of  performance  and  design  of  waterjets  was 
recently  presented  by  Kerwin  (2006). 

A boundary  element  method  (named  PROPCAV)  for 
the  analysis  of  cavitating  open  propellers  subject  to  a 
non-uniform  inflow  was  originally  developed  at  MIT 
by  Fine  (1992)  and  Kinnas  & Fine  (1992).  Since  then, 
the  method  has  been  improved  considerably  by  the 
Ocean  Engineering  Group  at  UT  Austin,  to  include 
mid-chord  back  and/or  face  cavitation,  modeling  of 
super-cavitating  and  surface-piercing  propellers  by 
Young  (2002)  and  Young  & Kinnas  (2001,  2003, 
2004);  unsteady  wake  alignment  and  developed  tip 
vortex  cavitation  by  Lee  (2002)  and  Lee  & Kinnas 
(2004,  2005);  modelling  of  cavitating  propellers  inside 
tunnels  and  of  cavitating  ducted  propellers  by  Lee  and 
Kinnas  (2006);  application  to  propeller  induced  mdder 
cavitation  by  Kinnas  et  al  (2007);  and  more  recently 
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coupling  with  an  integral  boundary  layer  solver  by  Sun 
and  Kinnas  (2006). 

In  this  work  we  present  an  extension  of  our  previous 
work  to  predict  the  performance  of  waterjet  propulsors, 
including  the  presence  of  sheet  cavitation  on  the  rotor 
and  stator  blades.  We  will  only  address  axial  flow 
pumps  subject  to  uniform  upstream  inflow  at  this  stage. 
The  interaction  between  rotor  and  stator  will  be  time- 
averaged  in  a similar  way  as  presented  in  Taylor  et  al 
(1998),  Kinnas  et  al  (2002),  Kerwin  et  al  (2006),  and 
Brewton  et  al  (2006). 


where  q(x,y,z)  is  the  total  velocity,  and  </>{x,y,z)  is 
the  perturbation  potential. 

Governing  equations  and  boundary  conditions 

Integral  equation  for  both  rotor  and  stator:  the 
perturbation  potential,  <j>(x,y,z) , at  any  point  p(x,y,z) 
located  either  on  the  wetted  rotor  or  stator  blades  and 
the  hub  surface,  SK  u Ss  cj  S„  , and  the  casing  surface. 
Sc , or  on  the  cavity  surfaces  of  the  rotor  or  stator, 
Src  Ssc  , has  to  satisfy  Green’s  third  identity. 


Figure  1:  Rotor-  and  Stator-Fixed  Coordinate  systems  and 
paneled  geometry  of  wateijet  components. 

FORMULATION 


A waterjet  geometry  with  the  related  coordinate 
systems  is  depicted  in  Figure  1. 

Assumptions 

- Inflow  at  waterjet  inlet  is  uniform  ( Uin ),  and  defined 
in  the  ship  fixed  coordinate  system. 

- Waterjet  rotor  rotates  with  a constant  angular 
velocity,??  . 

- Inflow  velocity 

Vm  (.r,  y,  z)  = Uin  (x,  r,  0)  in  the  case  of  a ship  fixed 
coordinate  system 

Vjn  (.v,  y,  z)  = Uin  (.r,  r.  O'-  ml)  -mxx  in  the  case  of  a 
rotating  coordinate  system 

- Fluid  is  inviscid,  and  the  flow  is  irrotational  and 
incompressible. 

q(x,  y,  z ) = Vln  (,r,  y,  z)  + V <f>(x,  y,  z ) 


Rotor/Stator  Probiem 


A* 


Figure  2:  Schematic  showing:  (a)  the  combined  rotor/stator 
problem,  (b)  the  rotor  problem,  and  (c)  the  stator  problem. 
Problems  (b)  and  (c)  are  solved  in  sequence  with  the  (time- 
averaged)  effects  of  one  on  the  other  being  accounted  for  in 
an  iterative  sense. 
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Integral  equation  for  wateriet  rotor: 
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where  the  subscripts,  q and  p,  correspond  to  the 
variable  point  and  the  field  point,  respectively. 
G{p,q)  = 1/ R(p]q)  is  the  Green  function,  where 
R{p\q)  is  the  distance  between  the  field  point  p and 
the  variable  point  q . nq  is  the  unit  normal  vector 


pointing  into  the  flow  field.  A</>m  and  are  the 
potential  jumps  in  the  trailing  wake  sheets  shedding 
from  either  the  rotor  or  the  stator  trailing  edge, 
respectively. 
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Where  4mf>RS  (x)  are  the  circumferentially  averaged 
values  of  the  induced  potentials  on  rotor,  hub  and 
waterjet  casing  due  to  stator,  defined  as  follows: 


2k0rs(x)  = J 
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Integral  equation  for  wateriet  stator: 


In  the  above  equation,  the  potentials  should  also  be  a 
function  of  time,  since  the  interaction  between  rotor 
and  stator  is  unsteady  in  nature.  The  above  equation 
can  be  applied  with  respect  to  the  rotating  coordinate 
system,  and  in  that  case  the  stator  is  a moving  surface 
and  the  appropriate  kinematic  boundary  condition  must 
be  applied  on  it,  or  with  respect  to  the  ship-fixed 
coordinate  system,  and  in  that  case  the  rotor  is  a 
moving  surface.  It  should  be  noted  that  the  value  of 

— for  each  component  is  independent  of  the 

dn 

coordinate  system. 

In  the  present  approach  we  will  consider  the 
circumferentially  averaged  effect  of  each  device  on  the 
other.  In  addition  we  will  solve  the  rotor  with  respect 
to  the  rotating  coordinate  system  and  the  stator  with 
respect  to  the  fixed  coordinates  system.  The  complete 
rotor/stator,  the  rotor,  and  the  stator  problems  are 
depicted  in  Figure  2. 
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+ 4 n<f>SR  (x) 


Where  ^x<f>SR  (x)  are  the  circumferentially  averaged 
values  of  the  induced  potentials  on  the  stator,  hub  and 
waterjet  casing  due  to  the  rotor,  defined  as  follows: 


2x</>sr  (x)  = | 
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Boundary  conditions: 


1 . The  flow  on  the  wetted  parts  of  rotor,  stator,  hub, 
and  casing  surfaces  should  be  tangent  to  the  wetted 
surfaces. 
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-Vjn  ■ n = -Ujn  ■ n + (m  x x)  • n;  on  the  rotor 
d<f> 

dn 

-Vm  ■ n = -Uin  -n  ; on  the  stator 

2.  The  Kutta  condition  requires  that  the  fluid  velocities 
at  the  rotor  and  stator  trailing  edge  are  finite. 

| V <f>{x,  y,  z)  |<  oo  at  rotor  or  stator  trailing  edge 

An  iterative  pressure  Kutta  condition  is  implemented 
as  described  in  Kinnas  & Hsin  (1992). 

3.  The  cavity  closure  condition  implies  that  the  cavity 
has  to  be  closed  at  its  end.  Since  the  cavity  planform  is 
unknown,  the  boundary  value  problem  is  solved  at  the 
given  cavitation  number  by  using  a guessed  cavity 
planform  which  may  not  be  closed  if  the  pressures  on 
the  cavity  planform  do  not  correspond  to  the  given 
cavitation  number.  The  Newton-Raphson  iterative 
method  is  adopted  to  find  the  correct  cavity  extent 
which  satisfies  the  cavity  closure  condition  at  the  given 
cavitation  number  (Kinnas  & Fine,  1993). 

4.  The  dynamic  boundary  condition  on  the  cavity 
siuface  requires  that  the  pressiue  on  the  cavity  surface 
is  constant  and  equal  to  the  constant  pressure,  pc.  inside 
the  cavity.  By  manipulating  the  Bernoulli’s  equation  in 
the  rotating  or  the  fixed  coordinate  system  in  terms  of 
the  cavitation  number,  the  total  velocity,  q, , on  the 
cavity  surface  can  be  given  as  follows. 

- i2  , , dd> 

Uin  + a>  r - 2 gy  - 2—;  on  rotor 
1 dt 

. „ -d</> 

^J-2gys-2— ;on  stator 
dt 
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Pin  d + < 


Figure  3:  Local  coordinates  system  at  each  panel.  The  shown 
s and  v axes  are  tangent  to  the  blade  surface.  Axis  w is  also 
tangent  to  the  blade  surface,  and  perpendicular  to  the  s axis. 

The  -2 dtfildt  term  in  the  Bernoulli  equation  will  be 
zero  in  the  case  of  axi-symmetric  inflow.  Thus,  for 
simplicity,  this  term  will  be  omitted  in  the  next 
equations,  even  though  it  is  still  present  in  the  code 
(where  it  will  also  become  zero  since  the  inflow  is  axi- 
symmetric  and  d/dt  = 0 ). 


The  total  velocity  qt  can  be  expressed  in  terms  of  the 
directional  derivatives  of  the  perturbation  potential  and 
the  inflow  components  at  the  non-orthogonal 
coordinate  system: 


t h 
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where  r is  the  distance  from  the  axis  of  rotation;  g is 
the  gravitational  acceleration;  ys  is  the  vertical 
distance  from  the  horizontal  plane  through  the  axis  of 
rotation;  n and  D are  the  blade  rotational  frequency 
and  the  rotor  diameter,  respectively.  The  cavitation 
number,  <rn  , is  defined  as  follows: 


with  s and  v being  the  unit  vectors  corresponding 
to  the  coordinates  s (chordwise)  and  v (spanwise), 
respectively,  and  with  n being  the  unit  normal  vector 
to  the  cavity.  US,UV,  and  U„  are  the  components  of 

the  inflow  velocity  Vin  along  the  (s,v,n)  directions. 
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1/2 pn2D2  s 
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v lev  VJ 

for  stator 


Vn  = j + U„  j are  the  components  of  the  total 
velocity  vector,  qt , along  the  directions  (s,v,n) . 


where  Po  is  the  pressure  far  upstream  at  the  depth  of  Combining  the  equations  for  \q\  and  qt  , given 
the  shaft  axis,  and  p is  the  water  density.  q < 

previously,  — can  be  obtained  as  follows: 

ds 
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where  9 is  the  angle  between  s and  v , as  shown  in 
Figure  3.  Note  that:  cos 9 = s - v 

A Dirichlet  type  of  boundary  condition  on  <j>  is  derived 

by  integrating  the  equation  for  — : 

ds 


On  the  part  of  the  cavity  over  the  blade  surface: 
<j>{s,v)  = <^(0,v)  + 


and,  on  the  part  of  the  cavity  over  the  wake  surface  is: 
<j>+  ( s,u ) = </>{sTE,n) 


+r 

-[/  + . |a,|2  - 

2 

J sTE 

5 'WmM 

{ du  j 

The  potential  ^(0,v)  corresponds  to  the  potential  value 
at  the  cavity  leading  edge,  and  can  be  extrapolated  in 
terms  of  the  unknown  potentials  of  wetted  flow  panels 
in  front  of  the  cavity  detachment  location.  5 = sTE 
denotes  the  blade  trailing  edge.  The  variable  u in  the 
equation  for  <f>+(s,n)  corresponds  to  the  directional 
derivative  normal  to  (5,  n)  plane  on  wake  surface,  and 
the  superscript,  +,  represents  the  upper  side  of  the 
wake  sheet.  The  equation  for  the  potential  on  the  cavity 

includes  the  unknown  functions,  — and  those  terms 

dv 

are  determined  in  an  iterative  manner. 


5.  The  kinematic  boundary  condition  on  cavity  surface 
requires  that  the  substantial  derivative  of  the  cavity 
surface  has  to  vanish.  Once  the  boundary  value 
problem  is  solved,  the  kinematic  condition  is  utilized  to 
determine  the  position  of  the  cavity  surface. 

\4j  + Vt  •vl[«-/?O,v)]=0 


where  h(s,v)  is  the  cavity  thickness  normal  to  the 


blade  surface. 

By  substituting  the  gradient  in  terms  of  the  local 
directional  derivatives,  the  partial  differential  equation 
for  the  cavity  thickness  is  derived  as  follows: 

- sin  9 ■ K ]+  |j[Fv  - sin  9 ■ Vs  ] = Vn  cos2  9 

Where 


K =—+us 

ds 


K=^~+uv 

OV 


vr,=-r+u » 
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The  cavity  height  nonnal  to  the  blade  surface  can  be 
determined  by  solving  the  above  partial  differential 

equation  which  includes  the  solution  — of  the  integral 

dn 

equation  for  the  rotor  and  stator. 

The  cavity  height  on  wake  surface,  ( hw ),  when  the 
super-cavity  occurs,  is  similarly  determined  in  terms  of 
the  cavity  source  on  the  wake  surface: 

dK  1-  I df  d<t> 

ds  ' dn  dn 

6.  The  cavity  detachment  location  on  either  side  of  the 
blade  is  determined  iteratively  to  satisfy  the  smooth 
detachment  condition  as  described  in  Young  and 
Kinnas  (2001).  The  cavity  detachment  may  also  be 
determined  via  coupling  with  XFOIL  (Drela,  1989)  as 
described  in  Brewer  & Kinnas  (1997)  and  Sun  & 
Kinnas  (2006). 

7.  The  velocity  at  the  inlet  has  to  be  equal  to  that 
specified,  i.e.  q = Vin  , and  that  leads  to  the  Neumann 
type  of  condition  at  the  inlet: 


d(j> 

dn 


= 0 at  the  inlet 


Similarly,  at  the  outlet  we  should  have: 

= (?«*-?*.)•» 

If  Uont  is  the  axial  component  of  qout  we  get: 


d(j> 

dn 


d(j> 

dn 


■Ut-U 


at  the  outlet 


Assuming  that  Uout  is  uniform1,  we  can  determine  its 
value  from  applying  continuity: 

77  A -TT  .A  —s TJ  -TI 

^ in  , in  ^ out  ^out  ^ out  ^ in  , 

^ out 

where  and  Aout  are  the  casing  areas  at  the  inlet  and 
outlet. 


Interaction  between  rotor  and  stator: 


The  fluid  field  around  the  waterjet  propulsor  is  solved 
in  an  iterative  manner  by  solving  the  integral  equations 
for  rotor  and  stator  separately  and  by  considering  the 


1 As  shown  in  Choi  & Kinnas  (1998)  this  assumption  has 
negligible  effect  on  the  flow  around  the  propeller. 
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effect  of  each  component  on  the  other  by  using  the 
induced  potentials  from  one  to  the  other.  The  induced 
potentials  on  the  other  component  are  calculated  using 
the  integral  equation  for  4 mj>RS  (x)  and  4x0s//  (x) , and 
are  circumferentially  averaged  to  apply  on  the  control 
points.  The  problems  for  both  components  are  iterated 
until  the  forces  converge  within  a certain  criterion. 
Figure  4 shows  how  the  induced  potentials  due  to  the 
rotor  are  included  at  the  control  point  on  the  stator,  and 
the  averaged  potentials  on  the  stator  induced  by  each 
panel  on  the  key  blade  of  the  rotor  can  be  derived  as 
follows  (N  is  the  number  of  equally  spaced  elements 
over  an  angle  equal  to  the  angle  between  two  rotor 
blades): 

i=N 

2>' 

J _ i=l 


Finally,  it  should  be  noted  that  the  swirl  (tangential 
velocity)  induced  by  the  rotor  on  the  stator  (assumed  to 
be  post-swirl)  will  also  need  to  be  evaluated  (by 
averaging  circumferentially  the  tangential  velocity  with 
respect  to  the  angular  position)  and  then  included,  as  a 
velocity  term,  in  the  kinematic  boundary  condition  on 
the  stator  blades.  Thus  the  kinematic  boundary 
condition  on  the  stator  blades  must  be  adjusted  as 
follows: 

= -(U,„  +Mtan,57?)-h 

stator 


d</> 

dn 


Numerical  Implementation 

A constant  dipole  and  source  panel  method  (code  name 
PROPCAV)  is  utilized  in  order  to  solve  either  the  rotor 
or  the  stator  problem.  Representative  panel 
arrangements  on  the  various  components  of  a notional 
waterjet  pump  are  shown  in  Figures  5-7.  The  panels  for 
the  rotor  or  the  stator  problems  are  aligned  with  the 
geometry  at  the  tip  of  the  rotor  or  the  stator  blades, 
respectively,  as  shown  in  Figures  6 and  7.  Special  care 
has  to  be  taken  with  the  panel  arrangement  at  the 
junctions  of  the  blades  with  the  hub  or  the  casing,  as 
shown  in  Figure  8,  in  order  to  avoid  highly  distorted 
panels. 

It  should  be  noted  that,  at  this  stage,  we  assume  that  the 
tip  gap  of  the  rotor  is  equal  to  zero,  even  though  the 
actual  gap  is  usually  under  1%.  As  shown  in  Kerwin 
(2006),  where  the  orifice  equation  was  implemented  in 
the  case  of  a wing  close  to  a wall,  for  a gap  of  1%  the 
results  with  zero  gap  were  much  closer  to  those  from 
inviscid  theory  with  the  orifice  equation  implemented, 
than  to  those  from  inviscid  theory  where  the  actual  gap 
was  used. 

The  results  are  represented  in  a non-dimensional 
manner,  as  follows: 

J = the  advance  ratio 
n D 


where  uXsaSR  is  the  tangential  (swirl)  velocity  induced 

by  the  rotor  on  the  stator  control  points.  Please  note 
that  this  adjustment  is  not  required  for  the  kinematic 
boundary  condition  on  the  rotor,  since  the  stator  does 
not  induce  any  (circumferentially  averaged)  swirl 
upstream  of  it. 
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Where  T and  Q are  the  thrust  and  torque,  respectively, 
acting  on  the  rotor,  D and  n are  the  rotor  diameter  and 
rotational  frequency,  A the  potential  jump  at  the 
trailing  edge  of  the  rotor  or  stator  blade,  and  Urcf  is  the 
reference  velocity  defined  as: 

Uref=4um2+{0.1nnD)2 


Figure  4:  Schematic  showing  the  inclusion  of  the 
circumferentially  averaged  rotor  effects  on  the  stator  control 
points 
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In  the  case  of  cavitation  the  pressure  on  the  cavity 
should  be  equal  to  the  vapor  pressure  and,  according  to 
the  definition  of  the  cavitation  number  for  the  rotor,  the 
following  equation  should  be  valid  on  the  cavity: 

—Cp  — <rn 

In  the  case  of  pumps  it  is  also  customary  to  evaluate 
the  headrise  (rise  in  total  pressure  head)  from  the  inlet 
and  outlet  section.  The  authors  plan  to  evaluate  this 
headrise  within  the  context  of  their  method  in  the  very 
near  future. 


VERIFICATION  AND  VALIDATION 

Results  and  grid  dependence  studies  in  wetted  and 
cavitating  flow 

To  study  the  numerical  performance  of  the  current 
model  the  method  is  first  applied  on  a notional  waterjet 
pinup,  which  is  based  on  the  one  which  is  currently 
being  tested  at  Johns  Hopkins  University  with  support 
by  the  Office  of  Naval  Research. 

The  results  from  several  grid  dependence  studies  are 
presented  in  Figures  9-19.  Figures  9 and  10  show  the 
effect  of  the  stator  on  the  rotor  as  the  number  of 
iterations  between  the  two  increases.  Note  that  the 
iterative  process  converges  very  quickly,  within  the  1st 
iteration  in  this  case  (the  0th  iteration  corresponds  to  the 
rotor  solution  without  the  stator).  Figiue  1 1 shows  the 
effect  of  the  stator  on  the  predicted  thrust  and  torque 
on  the  rotor.  Figiue  12  shows  the  convergence  of  the 
predicted  circulation  distribution  on  the  rotor  with 
number  of  panels  on  the  rotor.  Figiues  13  and  14  show 
the  convergence  of  the  predicted  circulation 
distribution  (wetted  and  cavitating)  with  the  number  of 
panels  between  the  blades  in  the  circumferential 
direction.  Figiue  15  shows  the  convergence  of  the 
wetted  circulation  distribution  with  the  size  of  the 
increment  in  the  rotational  direction,  A0,  which  sets  the 
size  of  the  panels  in  the  axial  direction  on  the  trailing 
wake  of  the  rotor  as  well  as  on  the  hub  and  casing. 
Figures  16  and  17  show  the  wetted  and  cavitating 
pressure  distributions  at  different  sections  along  the 
span  of  the  rotor  blades.  Please  note  that  the  cavitating 
pressure  distributions  are  such  that  -Cp  = crn  on  the 
cavity  and  -Cp  < <j„  everywhere  else  on  the  blade. 
Figures  18  and  19  show  contour  plots  and  the  predicted 
cavities  on  both  sides  of  the  rotor  blades.  It  should  be 
noted  that  midchord  back  and  leading  edge  face 
cavities  are  predicted  in  this  case. 


Figure  5:  Paneled  geometry  of  a notional  waterjet  pump  with 
a 5-blade  rotor  and  a 7-blade  stator.  Viewed  from  upstream. 


Figure  6:  Paneled  geometiy  of  a notional  wateijet  pump  for 
the  rotor  problem  (the  trailing  wake  of  one  blade  is  also 
shown).  Viewed  from  downstream. 


Figure  7:  Paneled  geometry  of  a notional  waterjet  pump  for 
the  stator  problem  (the  trailing  wake  of  one  blade  is  also 
shown).  Viewed  from  downstream. 
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Figure  8:  Details  of  paneling:  on  the  casing  at  the  blade 
leading  edge  (top  left),  on  the  casing  at  the  blade  trailing 
edge  (top  right),  on  the  hub  at  the  blade  leading  edge  (bottom 
left),  and  on  the  hub  at  the  blade  trailing  edge  (bottom  right). 


Figure  9:  Convergence  of  wetted  circulation  distribution  on 
rotor  with  number  of  iterations  for  a notional  waterjet  pump 
at  J=0.6  (the  0th  iteration  corresponds  to  the  rotor  solution 
without  the  stator) 


No.  of  iterations 


Figure  10:  Convergence  of  rotor  (wetted)  thrust  and  torque 
with  number  of  iterations  between  rotor  and  stator  for  a 
notional  waterjet  pump  at  J=0.6  (the  0th  iteration  corresponds 
to  the  rotor  solution  without  the  stator). 


Js 

Figure  1 1 : Predicted  rotor  thrust  and  torque  coefficients  with 
and  without  including  the  influence  of  the  stator  for  the 
notional  waterjet  pump  at  <Tn  =1 .0 
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Figure  12:  Convergence  of  the  predicted  wetted  circulation 
distribution  with  number  of  panels  on  the  rotor  blade  for  the 
notional  waterjet  pump  at  J=0.6 


r/R 


Figure  13:  Convergence  of  the  predicted  wetted  circulation 
distribution  with  the  number  of  panels  in  the  circumferential 
direction  between  the  blades  for  the  notional  waterjet  pump 
at  J=0.6 
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Figure  14:  Convergence  of  the  predicted  cavitating 
circulation  distribution  with  the  number  of  panels  in  the 
circumferential  direction  between  the  blades  for  the  notional 
waterjet  pump  at  J=0.6  and  crn  =1.0 


Figure  16:  Predicted  wetted  pressure  distributions  on  the 
rotor  for  the  notional  waterjet  pump  at  J=0.6. 
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Figure  17:  Predicted  cavitating  pressure  distributions  on  the 
rotor  for  the  notional  waterjet  pump  at  J=0.6  and  <Tn  =1.0. 


Figure  15:  Convergence  of  the  predicted  wetted  circulation 
distribution  on  the  rotor  with  panel  size  in  the  wake,  AO  , for 
the  notional  waterjet  pump  at  J=0.6  and  <7n  =1.0 


Cavrfy  <m  back  srde  (1-0.6) 


Figure  1 8:  Contour  plots  of  predicted  pressures  and  cavity 
patterns  on  the  suction  side  (back)  of  the  rotor  blades; 
notional  waterjet  pump  at  J=0.6  and  <Jn  = 1.0  , view  from 
upstream. 
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Figure  19:  Contour  plots  of  predicted  pressure  distributions 

and  cavity  patterns  on  pressure  side  (face)  of  the  rotor  blades; 

, . , . » Figure  2 1 : Thrust  and  torque  break-down  due  to  cavitation, 

notional  wateriet  pump  at  J=0. 6 and  = l.U , view  trom  ,.  . , . , ,,  , , ..  , ,,, 

J r r n predicted  by  the  present  method  when  applied  on  the  P1452 

downstream.  propeller  inside  the  D15  duct,  and  measured  by  Dyne  (1973) 

(From  Lee  and  Kinnas,  2006) 


Validation  for  ducted  propellers 


plication  to  the  CCDoTT  wateriet  i 


Some  correlations  of  the  present  method  when  applied 
to  ducted  propellers,  are  presented  in  Figures  20  and  21. 
Please  note  the  appreciable  effects  of  viscosity  on  the 
predicted  propeller  torque  shown  in  Figure  20.  Figure 
21  shows  the  effect  of  cavitation  on  thrust  and  torque 
breakdown. 


Fig.20:  Comparison  of  the  measured  and  predicted  blade 
forces  with  measured  on  PI 452  propeller  and  D 15  duct 
(Dyne,  1973).  The  viscous  results  have  been  produced  via 
interactive  coupling  of  the  present  method  with  XFOIL  on 
the  propeller  blades  (from  Kinnas  et  al,  2007) 


The  method  was  then  applied  on  a waterjet  designed  by 
CCDoTT  (Center  for  the  Commercial  Deployment  of 
Transportation  Technologies)  of  California  State 
University,  Long  Beach,  and  developed  in  conjunction 
with  CDI  Marine.  A side  view  of  the  waterjet  is  shown 
in  Fig.  22  and  its  main  characteristics  are  given  in 
Table  1.  The  CCDoTT  waterjet  was  tested  in  the  24- 
inch  water  tunnel  at  the  David  Taylor  Model  Basin  at 
NSWCCD  as  reported  by  Brewton  et  al  (2006),  who 
also  made  calculations  using  the  commercial  RANS 
solver  Fluent. 

The  paneling  on  the  rotor  and  stator  blades,  the  hub 
and  the  casing  are  shown  in  Figure  23.  The  predicted 
pressures  on  the  rotor  and  stator  blades  at  design  flow 
rate  are  shown  in  Figures  24-26.  Finally,  the  predicted 
vs.  measured  torque  on  the  rotor  blades  is  shown  in 
Figure  27.  Note  that  the  predicted  torque  is  lower  than 
measured  even  at  the  design  flow  rate.  The  authors 
think  that  due  to  the  rounded  trailing  edge  of  the  rotor 
blades  the  effects  of  viscosity  might  be  significant,  but 
at  this  stage  the  effects  of  viscosity  on  the  torque  were 
evaluated  by  just  applying  a uniform  friction 
coefficient  on  the  blades.  In  addition  the  authors 
assumed  a zero  gap  at  the  rotor  tip  even  though  the 
actual  gap  is  about  0.5%  of  the  rotor  radius.  Thus,  the 
effects  of  the  viscous  gap  flow  and  of  the  associated  tip 
gap  vortex  on  the  solution  were  not  accounted  for 
properly. 
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Figure  22:  A side  view  of  the  CCDoTT  waterjet  in  its  shroud 
(from  Brewton  et  al,  2006) 


Table  1 : Specifications  for  CCDoTT  pump  geometry  (from 
Brewton  et  al,  2006) 


Figure  23:  The  paneling  on  the  CCDoTT  waterjet.  The 
trailing  wake  of  one  of  the  rotor  blades  is  also  shown. 


Figure  24:  Predicted  pressures  on  suction  side  (back)  of  rotor 
blades  for  CCDoTT  wateijet  at  design  flow  rate.  Viewed 
from  upstream. 


Figure  26:  Predicted  pressures  on  stator  blades  for  CCDoTT 
watejet  at  design  flow  rate.  Viewed  from  downstream. 


Figure  25  : Predicted  pressures  on  pressure  side  (face)  of 
rotor  blades  for  CCDoTT  waterjet  at  design  flow  rate. 
Viewed  from  downstream. 
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Casing  Diameter 

0.1905  m 

Rotor  Diameter 

0.1895  m 

Rotor  Blades 

5 

Stator  Blades 

8 

Impeller  Chord  Length 

0.1516  m 

Stator  Chord  Length 

0.07389  m 

Rotor  Hub  Diameters 

0.0572  in  (forward); 
0.109  m (rear) 

Stator  Hub  Diameters 

0.109  m (forward); 
0.0508  m (rear  shaft) 

Rotational  Speed 

1700  rpm 

Mass  Flow  at  Design 

164.681  kg/s 
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Figure  27:  Predicted  by  the  present  method  and  measured 
torque  for  the  CCDoTT  wateijet  versus  the  ratio  of  the  flow 
rate  to  the  design  flow-rate. 


CONCLUSIONS 

A boundary  element  method  was  extended  to  predict 
the  performance  of  waterjet  pumps  subject  to  uniform 
inflow,  including  the  effects  of  sheet  cavitation.  The 
interaction  between  the  rotor  and  the  stator  was 
evaluated  in  an  iterative  manner,  by  considering  the 
circumferentially  averaged  effects  of  one  on  the  other. 

Preliminary  predictions  of  the  torque  on  a previously 
tested  waterjet  were  foimd  to  be  lower  than  those 
measured.  The  discrepancy  between  the  predictions 
and  the  measurements  might  be  attributed  to  the  effects 
of  viscosity  on  the  blade  pressures  and  torque  in  the 
case  of  roimd  trailing  edge  sections,  as  well  as  to  the 
effects  of  the  viscous  gap  flow  on  the  rotor 
performance. 

The  authors  intend  to  include  the  effects  of  viscosity 
on  the  rotor  blades  of  waterjets  by  extending  the 
method  of  Sun  & Kinnas  (2006).  They  also  plan  on 
including  the  effects  of  the  viscous  tip  gap  on  waterjet 
rotors  by  applying  the  orifice  equation  approach 
(Kerwin  et  al,  1987),  as  was  recently  implemented  in 
the  case  of  ducted  propellers  by  Kinnas  et  al  (2004),  as 
well  as  by  tracking  the  tip  gap  vortex  and  its  effect  on 
the  blade  loading  by  applying  the  method  presented  in 
Gu  (2006).  Systematic  correlations  with  the  detailed 
flow  field  velocity  and  pressure  measurements,  to  be 
reported  from  the  waterjet  tests  which  are  being 
performed  at  Johns  Hopkins  University,  will  also  be 
carried  out. 


Future  challenges  include  the  application  of  the  method 
to  cavitating  wateijets  (e.g.  correlation  of  predicted  to 
observed  cavity  planforms  on  rotor  or  stator  blades  and 
prediction  of  rotor  thrust  and  torque  breakdown  due  to 
cavitation),  and  the  prediction  of  the  unsteady 
cavitating  performance  of  wateijets.  The  latter  task  will 
allow  for  the  prediction  of  the  performance  of  waterjets 
which  are  often  subject  to  non-axisymmetric  inflows 
(due  to  the  wakes  of  the  waterjet  components 
upstream  of  the  rotor),  and  will  involve  the  evaluation 
of  the  3-D  (time-averaged  or  unsteady)  effective  wake 
to  the  rotor  or  stator  blades  via  coupling  with  a method 
(e.g.  a RANS  solver)  which  solves  for  the  global  flow 
inside  the  waterjet  including  at  least  the  inlet,  the  outlet, 
as  well  as  a part  of  the  hull  in  the  vicinity  of  the 
waterjet. 
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DISCUSSIONS  AND  AUTHOR’S  REPLIES 

NAME  OF  DISCUSSER: 

Henk  de  Koning  Gans 

QUESTION: 

The  wake  of  the  rotor  intersects  the  stator  so  a small 
slice  of  the  wake  is  inside  the  stator.  The  question  now 
is:  has  the  slice  of  the  wake  inside  the  stator  influenced 
the  boundary  condition  at  the  stator  collocation  point? 
And  has,  therefore,  the  inner  slice  of  wake  to  be 
subtracted? 

AUTHOR’S  REPLY: 

In  the  present  work  the  circumferentially  averaged  (or 
time  averaged)  effect  of  the  rotor  wake  on  the  stator  is 
considered.  Thus,  excluding  the  wake  elements  which 
are  inside  the  stator  blade,  should  have  a negligible 
effect  on  the  circumferentially  averaged  values.  In  the 
following  paper  by  Kinnas  et  al,  where  the  interaction 
between  a propeller  tip  vortex  and  a rudder  was 
addressed  in  an  unsteady  manner,  the  vortex  looped 
around  the  leading  edge  of  the  rudder,  and  was  not 
allowed  to  intersect  it. 

Kinnas,  S.A.,  Mohhamed,  F.,  and  Lee,  H.  “Propeller 
Tip  Vortex/Rudder  Interaction”,  Twenty-sixth 


Symposium  on  Naval  Hydrodynamics..  September 
2006,  Rome,  Italy. 
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NAME  OF  DISCUSSER: 

Ki-Han  Kim 

QUESTION: 

Where  did  you  apply  the  Kuttta  condition  for  the  blade 
with  a round  trailing  edge?  The  effect  of  the  location  of 
Kutta  condition  application  on  the  KT  and  KQ  for  a 
round  trailing  edge  would  be  significant. 

AUTHOR’S  REPLY: 

In  this  paper,  the  application  of  Kutta  condition  to  a 
round  trailing  edge  blade  is  not  addressed.  Instead,  the 
round  trailing  edge  is  modified  to  be  a sharp  one  by 
reducing  the  blade  thickness,  and  an  iterative  pressure 
Kutta  condition  is  applied  at  the  sharp  trailing  edge. 
The  authors  agree  with  Dr.  Kim  that  in  the  case  of  a 
round  trailing  edge  the  location  where  the  Kutta 
condition  is  applied  can  have  a significant  effect  on  the 
circulation  distribution,  and  thus  on  the  predicted  thrust 
and  torque.  However,  once  the  inviscid  solution  is 
coupled  with  the  boundary  layer  solver,  the  effect  of 
the  location  of  the  application  of  the  Kutta  condition 
on  the  final  viscous  solution  should  be  less  significant. 

-I 

NAME  OF  DISCUSSER: 

M.  Kinzel 

QUESTION: 

In  XFOIL,  do  you  allow  the  boundary  layers  to 
develop  over  the  cavity  or  do  you  allow  a slip  or  zero 
shear  stress  condition? 

AUTHOR’S  REPLY: 

For  sheet  cavitation,  the  water/vapor  interface  is 
treated  as  a constant  pressure  (free)  streamline.  The 
boundary  layers  are  allowed  to  develop  over  the  cavity 
surface  by  setting  zero  shear  stress  on  the  top  of  the 
cavity  surface,  i.e.  by  forcing  the  friction  coefficient  Cf 
to  be  0. 

A 

NAME  OF  DISCUSSER: 

J.  Young 

QUESTION: 

Is  the  presence  of  the  stator  accounted  for  in  the  wake 
alignment  with  the  rotor,  and  visa  versa? 

AUTHOR’S  REPLY: 

In  the  results  presented  in  this  paper,  there  was  no 
wake  alignment  performed,  and  the  rotor  and  stator 
wakes  were  set  at  the  same  pitch  as  the  corresponding 
blade  pitch. 
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